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ABSTRACT 

Chemical features of the local stellar disk have firmly established that long-term, 
continuous star formation has been accompanied by a steady rate of accretion of low- 
metallicity gas from the halo. We now argue that the recent discovery of an enhanced 
deuterium (D) fraction in the Galaxy is consistent with this picture. We consider two 
processes: the destruction of D in the interior of stars (astration) and the supply of 
nearly primordial D associated with the gas infall. Conventional Galactic chemical 
evolution models predict a monotonic decrease in D/H with time with a present-day 
D/H abundance which is much lower than the local value recently revealed. This 
predicted feature is the result of high levels of deuterium astration involved in the 
formation of the local metal-enhanced disk. Here we propose a new channel to explain 
the observed enhancement in D/H. Our model, which invokes ongoing gaseous infall 
and a star formation rate that declines over the past several Gyr, predicts that the D 
astration is suppressed over the same time interval. 
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1 INTRODUCTION 

During the past decade, there has been a steady improve- 
ment in observations of deuterium to hydrogen (D/H) in the 
universe. The measurement accuracy is now good enough to 
provide insight into Galaxy evolution. First of all, the pri- 
mordial D/H ratio (D/H) p is confined to a narrow range 
of ~(2.6-2.8)xl0" 5 . This range covers the analysis result 
of data from the Wilkinson Microwa ve Anisotropy P robe 
(WMAP) that gives 2.601° ^ xl0~ 5 dCoc et alj [2004 ) or 
2.75t° it xl0 ~ 5 ifCvburt et al1l2003l ), depending on the dif- 
ferent light-element nuc lear reaction rate ( also 2.57±0.15 
from the WMAP 5 data; iDunklev et aJj|200g ). and the mea- 
surements of distant m etal-poor damped Lyct systems that 
lead to 2.82jl°;f° x 10~ 5 (|Pettini et al.ll200sl . see also O'Meara 
et al. 2006). Recent Far- Ultraviolet Spectroscopic Explorer 
(FUSE) observations reveal that the present-day D /H abun- 
dance (D/H)o for interstellar matter (ISM) in the Galac- 
tic disk is surprisingly as high as ^ 2.31± 0.24xl 0~ 5 af- 
ter correcting for dust depletion (|Linskv et al.|[200^) . Some 
authors claim there are residual effects due to dust de- 
pletion whic h can account for the large scatter in th e ob- 
served D/H (|Steigman et al.ll2007l ; [Ellison et al.ll2007l ). But 
the high D/H abundance of 2.2±g"J x 10" 5 in the warm neu- 
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tral medium (WNM) of the lower Galactic hal o, which orig- 
inate s in the disk and is elevated into the halo ( Savage et al.l 
2007), strongly supports the claim bv lLinskv et al.l (|2006). 



The difference between (D/H) p and (D/H)o is defined 
as a deuterium astration factor /d=(D/H) p /(D/H)o. The 
value of fd should exceed unity because it is believed that 
all D is produced during Big Bang nucleosynthesis, and 
any subsequent destruction process of D occurs in the in- 
terior of stars. Thus, the value of fd is determined by the 
fraction of matter which has never been cycled through 
stars in the present ISM. Accordingly, fd is closely re- 
lated to the degree of heavy-element enrichment through 
star formation. Our expectation is that the D abundance 
in the ISM is anti-correlated with the metallicity such that 
the well-enriched (~solar metallicity) local disk demands a 
high fd. Indeed, conventional Galactic chemical evolution 
(GCE) models have been applied to th e evolution of D 
so far by several auth ors (|Prantzosl [l996 ; Tosi et al. I ll998l : 
iMatteucci et al.lll999r), and pre dict f d = 1.6 ± 0.2 for the 
local disk ( Romano et al.|[2006h . In other words, the local 
abundances of heavy elements such as Fe or a-elements in 
the ISM as well as long-lived stars imply fd ~ 1.4 at the 
least, that is equivalent to (D/H)o ~ 2.0xl0~ 5 . On the other 
hand, the high (D/H)o suggested by the recent observations 
implies a smaller fd. For the values of (D/H) p =2.75x 10~ 5 
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and (D/H)o=2.31xl0 -5 , we obtain / d =1.19, which is in- 
compatible with the theoretical prediction by GCE models. 

The D/H abundances which bridge between the primor- 
dial (D/H) p and the present-day (D/H)o provide crucial in- 
formation on the evolution of D. The D/H abundance in the 
protosolar cloud provides crucial information at a look-back 
time of ~ 4.6 Gyr. The D/H abundances deduced from the 
Jovian at mosphere rtLellouch et al.ll200j ) and from the so- 
lar wind l|Geiss fc Gloecklerll 19981 ) are estimated to be 2.1± 
0.4 (0.5) xlO -5 . This result is broadly consistent with the 
present-day local value. If these values are representative of 
the local disk, we are led to believe that D astration involved 
with star formation has turned off for the last several Gyr. 

In fact, as the process to change the amount of D in 
the ISM, the Galactic disk involves an enrichment process 
through an infall of low-metallicity gas from the halo, which 
brings (nearly) primordial D, besides its destruction through 
star formation. This infall is required by the GCE models 
to reproduce the ab undance dist ribution function (ADF) of 
local disk stars (e.g.. |Pagellll997t ). In addition, we cite high- 
velocity clo uds (HVCs) as the e vidence for the prese nt-day 
infall (e.g.. IWakker et al.|[l999l ; lLockman et allhoOSl ). The 
infall operates in the GCE models so as to enrich the ISM ef- 
ficiently in the early phase to suppress the number of metal- 
poor stars, and maintain the star formation activity by con- 
tinuously suppling the gas. This established view predicts a 
monotonic increase (decrease) in heavy elements (D) in the 
ISM up to the solar abundances ((D/H) £ 2.0xl0~ 5 ). In 
this way, the evolution of D is tightly bound to the produc- 
tion of metals in star formation. 

But what if a significant fraction of the heavy elements 
in th e disk does not originate from in situ star formation 
(e.g. lOppenheimer et alj |2010h ? This represents a radical 
departure from conventional GCE models and has recently 
been discussed in detail. Here large-scale winds from the 
Galactic bulge, which entrain a large amount of heavy ele- 
ments, enrich the disk. The hypothesis well expla ins the ori- 
gin o f metal-rich stars in the solar neighborhood (|Tsuiimotd 
120071 ) and the mechanism of the t ime-evolution of abundance 
gradient (|Tsuiimoto et al1l2010T ). In this scheme, a part of 
heavy elements trapped in long-lived stars as well as in the 
gas of the disk is not related to the local star formation but 
the activity of the bulge. 

In this paper, we shed light on the evolution of deu- 
terium in the framework that the Galaxy disk has evolved 
through "wind + infall", together with the renewed con- 
straint from the star formation history in the local disk. 
Then we show that the combined phenomena of winds that 
entrain heavy elements and infall that contains primordial 
D will make the Galaxy disk feature withstand the test of 
observations. 



2 WIND SCENARIO 

In this section, we summarize the proposed wind scenario, 
which is a key ingredient in the making of a hea vy-element 
featur e in t he disk, based on the d iscussions in ITsuiimotol 
< |2007f ) and iTsuiimoto et al.l ( |2010f ). together with the up- 
dated observational results. 



2.1 Signature of wind enrichment on the Galaxy 
disk 

1. Metal-rich stars in the solar neighborhood 

There exist metal-rich stars beyond the solar metal- 
licity, the fraction of which is roughly 20%. In fact, its 
presence as well as its elemental feature poses a puzzling 
problem. The ADF of solar neigh borhood disk stars has 
been studied by many autho r s (e.g., Wvse fc Gilmorel Il995l ; 
iRocha-Pinto fc Made] 1 19961 ; iNordstrom et all 120040 . and 
the location of its peak at [Fe/H] pea k~ -0.15 has been firmly 
established. On the other hand, the metal-rich end of the 
ADF extends at least to [Fe/H]~ +0.2. Spectroscopic ob- 
servations of elemental abundances for metal-rich disk stars 
|Feltzing fc Gustafssonlll998l ; iBensbv et alj|2005l ) have con- 
firmed that chemical enrichment in the s olar neighb o rhood 
continued unabated until [Fe/H]~ +0.4. ITsuiimotol (120071 ) 
has claimed that the simultaneous reproduction of both the 
presence of stars with [Fe/H] £ +0.2 and [Fe/H] poa k <0 is 
hard to realize through the conventional scheme of local en- 
richment via low-metallicity gaseous infall from the halo, 
and that the presence of supersolar stars is crucial evidence 
for enrichment by winds from the bulge. It should be of 
note that the seeming making of superso l ar stars by GCE 



models thus far (e.g. , Yoshii et 



supersolar sta rs by VjQil / 
al] 1 19961 ; lYin et all 120091 ; 



iKobavashi fc Nomotol [2009 ) is an end result of convolution 
of the ADF with a dispersion of 0.1 - 0.2 dex in [Fe/H]. 

The convincing signature of wind enrichment is seen in 
their elemental features. Nearby metal-rich stars are rep- 
resented by an upturn i n [q/Fe] ratios exc e pt for [O/Fe] 
again st increasing [Fe/H] i|Reddv et al.l 120031 ; IBensbv et all 
2005). Such upturns can also be seen in other elements, such 
as odd elements (Na, Al), and some iron-group elements (Ni, 
Zn). It is incorrect to attribute these upturns to metallicity- 
dependent SN II yields, especially because this feature is ob- 
served with both even- and odd-numberd elements. These 
features can be interpreted as a r eflection of the e lemen- 
tal feature of winds from the bulge (|Tsuiimotoll2007h . Since 
the [a/Fe] ratios start to upturn from [Fe/H]~ 0, the wind 
enrichment is expected to occur around 5 Gyr ago, corre- 
sponding to the time of sun's birth. From the above dis- 
cussion on elemental feature, the ga lactic fountain (e.g., 
IShapiro fc Fieldl 1 19761 : iBregmanl 11980) that may serve as a 
means for transporting metals to other locations in the disk 
is unlikely to contribute to acceleration of the chemical en- 
richment in the solar vicinity. 

On the other hand, radial mixing of disk stars due 
to resonant scattering with transient spiral waves pre- 
dicts that contaminants coming from the inner disk pop- 
ulate the solar neighborhood (Sellwood & Binney 2002; 
Roskar et al. 2008a,b; Schoenrich & Binney 2009; Sanchez- 
Blazquez et al. 2009), and allows the presence of metal-rich 
stars beyond the upper limit of a local chemical enrich- 
ment l|Roskar et al.|[2008bl ). However, stars orbiting near the 
Galactic center should exhibit a decreasing [ a/Fe] trend with 
increasing [Fe/H], as obser ved in the bulge jFulbright et al.l 
l2007tlMelendez et allboOSl ). Therefore, elemental features of 
nearby metal-rich stars seem at odds with the prediction of 
radial mixing as a mechanism of producing metal-rich stars. 
It should be of note that some authors claimed no upturning 
feature but a constant solar ratio for nearby metal-rich stars 
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l|Chen et al.ll2008h . However, none of the existing models in- 
cluding radial mixing supports this result. 
2. Time evolution of the ab undance gradient 

iTsuiimoto et alj |2010l ) has found that the radial [Fe/H] 
gradient roughly from J?gc=4 to i?GC = 14 kpc (-Rgc: Galac- 
tocentric distance) has flattened out in the last several Gyr 
with a change in a slope of ~-0.1 dex kpc -1 to ~ -0.05 dex 
kpc -1 , by comparing the steep relic gradient observed in the 
open clusters with that observed in t he Cepheids (see also 
IChen et al.l 120031 ; iMaciel et all 120061 ). Chemical evolution 
models predict contradictory tim e evolution of the me tal- 
licity gradien t, i.e., a steepening I Chiappin i et al.ll200ll) o r 
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a flattening |Boissier fe Prantzos! 1 1999 : iHou et al l 120001 ). 
Putting aside these contradictions, existing GCE models 
predict a monotonic increase in abundances for each region, 
and the predicted gradient change is sm all compared with 
the observations. ITsuiimoto et al.l (|2010h has claimed that 
large-scale winds from the Galactic bulge are the crucial fac- 
tor for its mechanism. Their proposed scenario is (i) winds 
once set up a steep abundance gradient (~ -0.1 dex kpc -1 ) 
owing to the enrichment by heavy elements entrained in the 
winds, and (ii) later evolution leads to a flattening of abun- 
dance gradient through chemical evolution under an accre- 
tion of a low-metal gas from the halo. 

On the other hand, new analysis of planetary nebulae 
supports a slight steepening of the [O/H] gradient with time 
(|Stanghellini fe Haywood 2010). This contradictory result is 
interesting because it is likely that winds from the bulge 
are not so metal-rich in O. The metal-rich bulge stars ex- 
hibit the low [O/Fe] below the solar ratio |Fulbright et al.l 
120071 ). Abundance of the microlensed stars recently an- 
alyzed confirmed [Q /Fe]~ -0.2 for metal-rich bulge stars 
(|Bensbv et al.ll2010T ). The reduced O yield of SN II expected 
from these observations is probably due to efficient mass 
loss for massive stars, which prevents He and C from syn- 
thesizi ng into O within the m etal-rich environment in the 
bulge (|Mc William et al. 1 120081 ). Thus, O is not a good in- 
dicator for the wind enrichment that is predicted to cause 
the time evolution of the abundance gradient in our scheme. 
In other words, the time evolution of the [O/H] gradient 
could be possibl y discussed in the conve ntional GCE scheme 
as predicted by IChiappini et all j2001). Alternatively, it is 
possible that some effect by radial mixing which causes 
a steepening of abundance g radient |Roskar et all [2008b; 
ISanchez-Blazquez et al.1 12009) is seen in the evolution of O 
gradient owing t o little influence of the w i nd en richment. 

In addition, IStanghellini fe Havwoodl |2010l ) compiled a 
new data set of ages, [Fe/H], and distances for open clusters. 
From their samples, we make two groups: one is composed 
of young clusters whose ages are less than 1 Gyr, and the 
other is of old ones whose ages range within 5 ± 2 Gyr. Here 
we restrict both of samples to the ones with -Rgc <14 kpc. 
As a result, the least-square fitting gives a slope of -0.05 
dex kpc -1 for young clusters, and -0.14 dex kpc -1 for old 
clusters, respectively. We see a significant flattening of the 
[Fe/H] gradient for new open cluster data. 
3. Other signatures 

In the Galaxy, there is now strong evidence 
for lar ge-scale winds across th e elec tr omagnetic spec- 
trum iBland-Hawthorn fc Cohenl 120031 ; iFox et al.1 120051 ; 
iKeenev et alll2006l : lEverett et al.ll2010l ). Furthermore, there 
is tantalizing direct evidence that dust and gas from the 



inner Galaxy have b een transported to the solar neighbour- 
hood l|Clavtonlll997l ). In standard chemical evolution mod- 
els, the isotopes 29 Si and 30 Si become more abundant than 
28 Si as the disk ages. But pre-solar grains from meteorites 
show evidence that the local isotopes formed in material 
that had experienced more nuclear synthesis than material 
that came after the Sun's birth. 



2.2 Driver of winds 

A driver of large-scale winds would be a starburst in the 
Galactic bulge. Therefore, our scenario proposes that the 
Galactic bulge has a complex history involving at least one 
starburst episode around 5 Gyr ago in addition to the ini- 
tial burst that gave rise to the old population > 10 Gyr ago. 
To date, the star formation history of the Galactic bulge 
still remains an open question because we do not have the 
precise color- magnitude diagram (CMD) of the bulge ow- 
ing to the diffic ulty of discrimination of bulge stars from 
disk stars (e.g.. Irloltzman et al.l Il993: Fcl tzing fc Gilmord 
2000). There is go od evidence for recent star formation in 
the nuclear region l|Launhardt et alj2002h . In addition, from 
the perspective of bar formation, one might expect some 
star formation activity in the region of what is now the 
bulge, about 3 Gyr after the inner disk started to form - 
this would be at the time when the bar buckling occurs 
|Athanassoula fc Martinez- Valpuestall2009l ) - it would really 
stir up any remaining gas in the inner disk. 

It is difficult to conjecture how massive a starburst is 
likely to be from a viewpoint of a wind enrichment on the 
disk since there are too many uncertain factors to predict 
it. Here we take a simplified case and show that a starburst 
which produce ~ 5% of the bulge stars might be one candi- 
date, by checking supply and demand of Fe. For simplicity, 
the disk is here considered to be one zone represented by the 
solar circle. Let us take the masses of a bulge and a disk and 
the gas fraction of a disk as 2 x 10 10 M B , 4 x 10 10 Mp), and 
0.1. The IMF with a slope of x= -1.05 jMatteucci fc Brocatol 
1990; iBallero et al.ll2007l ) together with the star formation 
producing the fraction of 5% of bulge stars yields ~ 3 x 10 7 
SNe II that eject ~ 4 x 10 6 M Q of Fe, some fraction of which 
will be, though, retained within the bulge. On the other 
hand, the Fe abundance of gas in the disk should increase 
from [Fe/H]=0 to [Fe/H]=+0.3 by the wind enrichment plus 
the local enrichment by star fromation. It means that the 
winds should entrain roughly the same amount of Fe as that 
already present in the ISM to double the Fe abundance if we 
neglect the local Fe production. In the end, the total Fe mass 
necessary for the additional enrichment on the disk amounts 
to 4 x 10 10 x 0.1 x 10~ 3 (solar Fe abundance) ~4x 10 6 M Q , 
equivalent to the former estimate. If we assume that the 
wind enrichment is restricted to the inner disk including 
the solar vicinity together with a partial contribution to the 
winds from type la supernovae (SNe la) and hypernovae 
that eject high Fe masses, the necessary fraction of stars 
produced in a burst will significantl y decrease like < 1%. 

While numer ous stu dies (lOrtolani et all 1 19951 : 
Feltzing fc Gilmorel l2000l: fKuiiken fc Rich! 120021 : 



Zoccali et al.1 120031 ; IClarkson et al.1 l2008fi agree that 



the bulk of the bulge formed >10 Gyr, a minority popula- 
tion of stars as young as 5 Gyr cannot be ruled out, and a 
minor burst of star formation as recent as 5 Gyr may have 
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taken place in the bulge. The fraction of stars produced 
by such a minor burst must be small like ~ 5% as implied 
by an above simplified estimate so as not to contradict 
the present CMDs, but in practice, how much of stars the 
starburst was likely to produce remains quite unclear. As 
well, the age of a burst claimed here is not definitive due to 
a lack of determinant factor of its age-dating. The results 
from microlensed dwarf and subgiant stars in the Galactic 
bulge showing the averag e age of ~ 7 Gyr fo r metal-rich 
([Fe/H]> 0) bulge stars l|Bensbv et al.l l201fj ) may be in 
favor of our hypothesis. 

Recently, some authors claim that momentum- 
driven winds are an efficient mech anism for a power- 
ful wind from the ce n tral s ta rburst ([Murray et al.l 20051 ; 
iQppenheimer fc Pavel 120061 ; lOppenheimer et al l I20T0F . 
These strong winds could realize a widespread enrichment 
over the disk, possibly including the outer disk region, 
through hydrodyna mical interac t ion b e tween the winds and 
the gaseous halo (|Bekki et all [2009; Oppenheimer et al.l 
2010). In N-body hydrodynamic simulations, only high- 
resolution models can create inhomogeneous ISM sufficient 
for t he right treatment of win ds from star-forming region 
fe.g.. [Ceverino fc Kly pin 2009). Ongoing advanced numeri- 
cal calculations will certainly reveal a full view of the wind 
enrichment. 



2.3 Influence of an initial starburst in the 
Galactic bulge 

The main body of the Galactic bulge was considered to be 
formed in a starburst > 10 Gyr ago. It is natural that we 
expect the action of ravaging winds associated with this vast 
star formati on. Since the age of thin disk is as low a s ^ 
10 Gyr (e.g., Hansen et alj|2002t Idel Peloso et al.ll2005l ), the 
signature of initial wind enrichment is likely to be unseen 
in the thin disk stars. Instead, some fossil record should 
be imprinted in the formation of thick disk prior to thin 
disk. Indeed, recent study claims that ages of the thick disk 
stars are >8-9 Gyr with a span of ~3 Gyr l|Bensbv et al.l 
120071 ), which suggests a similar formation epoch between 
thick disk and bulge. Here we discuss compelling evidence for 
the contamination of a gas used for the thick disk formation 
by heavy elements entrained in the winds. 

The ADF of thick disk shows a sharp increase from 
[Fe/H]~ -1.3 to the peak l ocated at [Fe/H]~ -0.7 . This 
feature is first reported by IWvse fc Gilmorg fl995) with 
a small number of but complete samples, and the loca- 
tion of its peak is confirmed by a huge SPS S data base 
of thick disk stars (|Allende Prieto et al.l l2006h . These ob- 
servational results suggest that a material for the proto- 
thick disk was pre-enriched up to [Fe/H]~ -1.3 at > 10 
Gyr ago. Our claim is that the mechanism of pre-enrichment 
is attributable to the wind occurrence triggered by an ini- 
tial starburst in the Galactic bulge. Petailed results of not 
only the APF but also the evolution of [a/Fe] calculated 
with the wind models will be presented in Tsujimoto (2010, 
in preparation). In addition, similarity of stellar elemen- 
tal features between bulge and thick disk r ecently revealed 
[[Melendez et al.1 |200S| ; IBensbv et ail I2010T I seems to have 
some connections to the wind phenomena during their over- 
lapping formation episode. 

From this viewpoint, metal- weak thick disk stars can be 



discussed. IChiba fc Beers! ([20001 ) found that some fraction 
of low-metallicity stars in the range of —2.2 ^[Fe /H]^ — 1 
belong to the thick disk population. According to ICaimmil 
(2008 ), addition of the se stars to the APF deduced by 



IWvse fc Gilmord |l995) does not form a metal-poor tail. 
Rather, it produces another distinct distribution for [Fe/H]< 
— 1.4. In our view, metal-weak thick disk stars are regarded 
as the stars that are formed prior to the significant influence 
of metal-rich winds from the bulge on the proto-thick disk. 



3 MODEL AND OBSERVATIONAL 
CONSTRAINTS 

3.1 "wind+infall" model 

We incorporate the enrichment process by large-scale winds 
into the standard GCE model. The basic picture is that the 
disk was formed through a continuous low-metal infall of ma- 
terial from outside the disk region, and for a specific period 
heavy elements carried by winds from the bulge dropped and 
enriched the disk. The period of the occurrence of winds is 
set to be 4-6 Gyr ago so as to reproduce the chemical fea- 
tures of metal-rich solar neighborhood stars as well as the 
obser ved flattening of abund ance gradient in the last several 
Gyr ([Tsujimoto et alj|20ich . The abundances of the winds 
are set to be identical to the ejecta of type II SNe (SNe II) 
such as [Fe/H]=+0.4, [Mg/Fe]=+0.4, and [O/Fe]=-0.2. A 
low O abundance in the wind is implied by the bulge abun- 
danc e j|Fulbright et al.ll2007l ; lTsuiimotdl2007l ; IBensbv et all 
|2010|) . Petails of the mo dels including th e equa tions to be 
solved are described in iTsuiimoto et al.l ([20105). The key 
model ingredients, infall rate and star formation rate (SFR) , 
are presented in the following separate subsections. On other 
model inputs, their individual descriptions are simply men- 
tioned below. 

The initial mass function (IMF) should be chosen care- 
fully in this study since it determines the returned frac- 
tion R, i.e., the fraction of stellar mass that is returned 
to th e ISM, which crucially infl uences the evolution of P 
(e.g., iProdanovic fc Fieldsl I2008T) . We adopt a power-law 
mass spectrum with a slope of -1.35, together with the mass 
range [0.05 Mq, 50 M@]. This combination is deduced from 
the theoretical analysis considering two conditions given by 
the chemical abundances of halo stars a nd the mass-to- 
lumin osity ratio in the solar neighborhood l|Tsuiimoto et all 
Il997h . Then, this adopted IMF is combined with the nucle- 
osyn thesis yields of SNe II a nd type la SNe (SNe la) taken 
from ITsuiimoto et al.l l|l995l ). For the occurrence frequency 
of SNe la, we assume that the fraction of the stars that even- 
tually produce SNe la for 3 — 8Mq is 0.05 with the l ifetime 
of SN la progenitors of 0.5-3 Gyr (lYoshii et al.lll996l ). Note 
that the returned fraction appears in the equations to solve 
the evolutions of gas and abundances (O, Mg, Fe) with no 
instantaneous recycling approximation. 

The evolution of P is described by a simple equation 
that is governed by the time changes of a gas fraction f g , a 
SFR ip(t), and an infall rate A(t); 



d(Df 9 ) 
dt 



= -DiP(t) + D p A(t) 



(1) 



where P in the winds is neglected since the winds mostly 
consists of the ejecta of SNe II with no P. The P abundance 
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Figure 1. upper panel: The model infall rate as a function of look- 
back time, together with th e observed presen t-day rate which is 
estimated from Complex C jThom et al ]|2008l . ; see, however, the 
discussion in the text). The model infall rate is normalized so that 
its integration over the whole age gives a total surface density of 
42 Mq pc - 2 . middle panel: The calculated star formation histo- 
ries. Two model cases are considered. One (solid line) involves a 
setup of a decreasing SFR for 0< t <4 Gyr so as to b e compatible 
with t he observed tendency (dashed line) derived bv lFuchs et al.l 
(2009). The other (dotted line) is the one with an assumption 
of a constant SFR, coeficient. Both models are normalized to be- 
come the present-day stellar surface density of 35.5 Mq pc -2 . 
The observed SFR at the present time is ded uced from the mas s 
growth of the stellar disk in the last 1 Gyr (see lFuchs et al . 2009). 
lower panel: The evolution of a gas fraction calculated with the 
decreasing SFR model, together with the present-day value. 



in the infall is assumed to be primordial since the metallic- 
ity ([Fe/H]) of the infalling gas is expected to be much lower 
than the present-day infall such as Complex C whose metal- 
licity is subsolar and D/H is estimated to be 2.2±0.7xl0 -5 
ijSembach et al.|[2004h . In the following, we will describe the 
key ingredients that determine the D/H evolution, together 
with the observational constraints to each of them. 



3.2 infall rate 

For the infall rate, we adopt the formula that is proportional 
to texp(— i/ti n ) with a timescale of infall ri n . The value of 
ri n is assumed to be 5 Gyr so as to rep roduce the observe d 
feature of metal-poor tail of the ADF (|Yoshii et al.lll996h . 
To compare with the observations, we give an infall rate in 
units of Mq yr -1 kpc -2 , that is shown on the upper panel 
of Figure 1. Here we normalize this function so that its in- 
tegration over 14 Gyr sums up to today's total surface den- 
sity in the local disk. Today's density is estimated to be 42 
Mq pc -2 from the surface density of the stellar disk (35.5 



Mq pc~ 2 ; lGilmore et al.|[l995l ; iFlvnn et al 



2006) and of the 



interstellar gas (6.5 Mq pc - ; [DameJ [1993J) . This together 
with our infall function predicts the present rate of 0.00185 
yr -1 kpc -2 . On the other hand, Complex C giv es an associ- 
ated mass flow of 0.00067 - 0.003 yr -1 kpc -2 l|Thom et all 
2008). In the figure, we adopt these values as an observation- 
ally implied range. However, the average rate given by the 
entire HVC system may be higher by a factor of ~ 2 as pre - 
dicted so far (e.g.. IChiappini et al.ll200ll ; iPeek et al.ll2008f ). 
Still, our prediction for the present-day rate falls within the 
observationally possible range. Furthermore, it should be 
of note that a significant number of ionized high-velocity 
gas detected in the vicinity of the Gal axy must be a factor 
to increase the present-day rate (e.g., Scmbach et al. I l2003l ; 
iFox et al.|[200l ). 



3.3 star formation rate 

The SFR is assumed to be proportional to the gas fraction 
f g with a coefficient v (Gyr -1 ). In the standard GCE 
models, the value of v is set to be constant over the age of 
a galaxy. On the other hand, recent several works reveal 
the star formation hi story in the local disk especially for 
the l as t several Gyr ([Hernandez et alj |2000|; IVergelv et al.l 



2002; Ide la Fucntc Marcos & dc la F uente Marcos! 12004 



ICignoni et al. 2006; Fu chs et al.ll2009l ). The common feature 



claimed by them is that the SFR has been clearly declining 
with time for the last several Gyr. For the purpose of its 
reproduction, we introduce the time-dependent v for the 
recent SFR so as to obey t he observed SFR- t relation. 
Using the relation found by iFuchs et alj l|2009l ) giving a 
decreasing rate of ~ -0.8 Mq Gyr -1 , we assume ^=0.6 for 
t < 10 Gyr and v=0.6x exp(-0.5 x (t - 10)) for 10< t < 
14 Gyr. The middle panel of Figure 1 shows the history of 
SFR where its normalization is done so that the present-day 
stellar density beomes 35.5 Mq pc -2 . For reference, the 
case with a constant v (=0.6) over the full age is indicated 
by dotted line, demonstrating that a constant v model 
cannot reproduce the observed SFR over the last Gyr. 

The present SFR is deduced from the mass growth 
of stellar dens it y ove r the last 1 Gyr. From the result 
of IFuchs et al.l 1 20091 ) with an associated uncertainty, it 
is estimated to be 0.48 - 1.1 Mq pc -2 Gyr -1 , which 
covers the results obtai ned from other studies including 
iRocha-Pinto et ail (2000). On the other hand, the predicted 
SFR is 0.82 Mq pc -2 Gyr -1 . We see a good agreement with 
the observations. 
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Figure 2. Features of chemical evolution in the solar neighborhood, left panel: Predicted abundance distribution function of thin disk 
stars against the Fe abundance. A good fit to the observation is realized by the model with winds (solid line) while the result without 
winds is shown by the dotted line. The calculated distributions are convolved using the Gaussian with a di spersion of 0.05 dex as t his is 
the minimum measurem ent error expected in the data. Open circles and crosses repr esent data taken from [ Edvardsson et al. I l ll993h and 
IWvse fe Gilmorel dl995T) . respectively. The model di stribution and the observ ed one bv lWvse fe Gilmorel lll995T ) are normalized to coincide 
with the total number of the sample stars used by [Edvardsso n et al. | il99$) .riaht panel: Correlation of [Mg/Fe] and [ O/Fe] with [Fe/H] 
for th in dis k stars with our win d model and no wind model superimposed. The open circles and dots are taken from [Edvardsson et al] 
il993l ') and lBensbv et alj {200 J) , respectively. 



3.4 gas fraction 

Finally, we see the evolution of gas fraction f g as a result of 
adopted infall rate and SFR since f g is a determinant factor 
to control the D/H abundance (the lower panel of Fig. 1). 
The predicted f g is compared with the observed value at the 
present time. An error in the observed f g mainly originates 
from a large uncertainty in the surface densit y of the ISM . 
Here we adopt its range from 6.5 Mq pc~ 2 (|Damel Il993h 
to 13 M© pc~ 2 (|Kulkarni fe Heiles]|l987l ). together with the 
stellar density of 35.5 Mq pc -2 . A gradual increase in f g 
over the last 4 Gyr while the SFR decreases seems at odds 
with the Schmidt-Kennicut law. However, a span in the SFR 
against a small range of f g in our result is compatible with 
the intrinsic scatter in the observed rel ation between th e 
SFR and the surface density of the ISM (|Kennicutt Ill99ct ). 
What primarily induces this implied low star formation ef- 
ficiency (SFE) irrespective of the steady supply of H I gas 
remains unclear, including the possible association with the 
low SFE deviated from the standard law in a low gas den- 
sity environme nt such as 5-10 Mp pc~ 2 o bserved for faint 
dwarf galaxies l|Rovchowdhurv et alj|2009l) a nd low surface 
brightness (LSB) galaxies (|Wvder et alj2009h . Interestingly, 
LSB galaxies exhibit about one-order change in SFR while 
those have similar gas densities. In addition, the Schmidt- 
Kennicutt law is valid only in a large scale such as > 1 kpc 
jOnodera et al.ll2010h . whereas the local SFR observation- 
ally covers tens of pc. 



4 HEAVY-ELEMENT FEATURES 

We proceed to the next issue on chemical features of the 
local disk. Recall that we need to discuss the D/H evolution 
in the framework that gives a fully consistent explanation 
to the observed quantities on heavy elements such as Fe and 
a-elements. The left panel of Figure 2 shows the successful 



reproduction of the ADF by our model. A stall of chemical 
enrichment owing to the decreasing SFR over the last Gyrs 
is compensated by the effective wind enrichment. In other 
words, the activity of local star formation is not sufficient to 
provide all of Fe contained in the existing long-lived stars. 
This is shown by the model case without the winds (dashed- 
line), which is skewed to the lower- met allicty, compared with 
the observations. 

This wind enrichment imprints its unique record on the 
elemental ratios of metal-rich stars. The winds are expected 
to exhibit the high [Mg/Fe] and low [O/Fe] ratios as men- 
tioned in the model description. Accordingly, as a result of 
their reflections, our model predicts an upturning feature 
of [Mg/Fe] for metal-rich stars through a rapid increase- 
gradual decrease loop in the [Mg/Fe]- [Fe/H] correlation, 
which is in good agreement with the observed data (the right 
upper panel). The different [O/Fc] evolution can be consid- 
ered in the same theoretical scheme with an assumption of 
a low [O/Fe] in the winds (the right lower panel). 



5 EVOLUTION OF D 

Now we are all set to start on the discussion of D/H evolu- 
tion. The predicted time-evolution of D/H is shown in Fig- 
ure 3. Here the pri mordial (D/H) p is set to be 2.75x10 
|Cvburt et al.| [2003). From its initial value, the D/H abun- 
dance starts to decrease with time, and subsequently its 
change turns to a gradual increase 4 Gyr ago. This turnover 
results from the supply of primordial D by an infall that 
overwhelms the destruction of D through star formation 
owing to its low rate. Then, it finally enables to achieve 
a present-day high D/H abundan ce of 2.3xl0~ 5 , ful ly con- 
sistent with the observed value |Linskv et al.l 120061 ). after 
passing D/H ~ 2.1 x 10~ 5 equivalent to the pro tosolar value 
l|Geiss fe Gloecklerl 1 19981 ; iLellouch et ail l200lh around 4.6 
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Figure 3. The evolution of D/H abundance in the local disk pre- 
dicted by the wind+decreasing SFR model (solid line). For ref- 
erence, the model case with a constant SFR coefficient is shown 
by the da shed line. Three ob served points are for the primor- 
dial value llCvbu rt et a l.ll2003h. the value for t he protosolar cloud 
llGeiss fc Gloecklerl 1998l; lLellouch et alfcOOll ). and the local ISM 
value l|Linskv et al.ll2006l ' . 



Gyr ago. For reference, the case with a constant v (=0.6) 
over the full age is indicated by dashed line. 

Our safe conclusion is that the D/H abundance has lit- 
tle evolution during the last several Gyr, as if the D as- 
tration through star formation has turned off. It's because 
the metallicity of an infall would change in acc ordance with 
the c osmic evolution of damped Lya systems l|Wolfe et al.l 
2005), and eventually become subsolar with a smaller D tha n 
the primordial value like Complex C (|Sembach et al.|[2004h . 
This view will lead to a flattening of D/H evolution rather 
than its gradual increase. 



6 UNCERTAINTIES IN OBSERVATIONAL 
CONSTRAINTS 

In this section, we briefly review the uncertainties in two 
key observational constraints that largely influence the re- 
sults presented here. As discussed below, neither of two has 
reached a firm consensus yet, thus further studies to better 
understand these issues, observationally and theoretically, 
should be still awaited. 



isochrones, lAumer fc Binnevl l|2009l ) obtain the deceasing 
SFR from the past to the present but at a very small rate. 

With no doubt, major progress on this issue will come 
from a huge database brought by GAIA in the next decade. 



6.2 present-day D/H abundance 

The controversial situation on the r eprese ntative value of 
(D/H)o is nicely reviewed in iTosil (|2010l . and references 
therein). Although most claims seem inclined to be above 
(D/H)o ~ 2 x 10~ 5 , that is beyond the predictions by 
the conventional GCE models, very low values such as 
(D/H) ~1.5xl0" 5 of the local bubble value or ~0.98x 10~ 5 
of the highest column density regions (|Hebrard et al.ll2005l ) 
can not be still discarded. The difficulty lies in the interpre- 
tation of a widespread variation in the observed gas-phase 
D/H ratios among local ISM regions. It seems reasonable to 
claim that the ISM exhibiting a high D/H such as ~2.2- 
2.3xl0~ 5 is a result of incompletely mix ed infall, i.e., a 
stron g influence by a low-metal infall (e.g.. Isteigman et al.l 
l2007f) since Complex C has a similar D/H of 2.2xl0" 5 . 
However, it is stemmed by a counterargument that their 
O a bundance is no t so low as that expected in an infalling 
gas l|Hebrardll2O10l ). 

Recently, analysis of Linsky's data by a Bayesian 
approach gives the be st estimate of (D/H)o ~2.0xl0 
jProdanovic et aT1l2010l ). A reduction of ~ 0.3 x 10 -5 draws 
the observed (D/H)o close to the compatibility with some 
of GCE models. However, considering an uncertainty in the 
adopted probability distributio n for depletion as w ell as an 
indication of high D in WNM jSavage et al.ll2007r i. further 
studies with more data is still necessary for the settlement of 
a true (D/H)o. Besides, recall that the observed D/H abun- 
dance only gives a lower limit of (D/H)o. 

T o date, a high (D/H)o as proposed by iLinskv et al.l 
(2006) is viewed with a skeptical eye owing to an incompat- 
ibility with the theoretical GCE scheme. Here we propose a 
new path to end up with (D/H)o ~2.3xl0 . Accordingly, 
we claim a need to have a fresh look at those observed high 
D/H, bearing in mind that they are within the predictions 
of a renewed GCE picture. 



6.1 recent star formation rate 



7 ROLES OF SFR AND WINDS 



In our models, we adopt a sharply decreasing SFR 
which cha nges a facto r of ~ 4 in the last 4 Gyr de- 
duced by iFuchs et al.1 (|2009p . Their result is broadly 
consistent with other several works using differ- 
ent methods dHernandez et al.l 2000l; Vergelv et al.l 



20021; Idela Fuente Marcos fc de la Fuente Marcos! 1200*41 
Cignoni et al. 20061 ) . However, at least t wo studies present 
a diffe rent picture for the recent SFR. iRocha-Pinto et aU 
(2000) derive the star formation history, using a chromo- 
sphcric age distribution of dwarf stars, and their major 
finding is that the disk has experienced a few bursts in the 
past, in which one of enhanced episodes of star formation 
occurred during the last 1 Gyr. On the other hand, from a 
new reduction of Hipparcos data combined with updated 



Finally, let's review each role of the sharply decreasing SFR 
and the metal-rich winds that are the main drivers of the 
successful results through the comparison of two models. 
Left panels of Figure 4 show the time evolution of [Fe/H], 
[Mg/Fe], and D/H calculated with the model in which the 
sharply decreasing SFR for the last 4 Gyr is incorporated 
into the GCE scheme. These witness that an inactive SF 
period makes the D/H evolution compatible with the obser- 
vation, while leaving the problem of an incomplete chemical 
enrichment. Then, the satisfactory enrichment is achieved 
by an introduction of the metal-rich winds into the model 
with the D/H feature unchanged (right panels). For refer- 
ence, in both models, the results calculated with a constant 
SFR coefficient are indicated by dashed lines. 
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Figure 4. Time evolution of [Fe/H], [Mg/Fe], and D/H. left pan- 
els: Solid and dashed lines represent the results by the model with 
the decreasing SFR and with a constant SFR coefficient, respec- 
tively, under the assumption of no winds, right panels: Same as 
in left panels but for an introduction of the winds. 

8 CONCLUSIONS 

As long as we stick to the conventional scheme of chem- 
ical evolution, recent observational implications of (i) the 
present-day high D/H abundance and (ii) little evolution of 
its abundance during the last several Gyr in the local disk 
are hard to believe. Accordingly, the discussion on it so far 
has been inclined to favor a much lower present-day D/H 
abundance in terms of an uncertainty in its estimate mainly 
due to dust correction. 

Here we present a new view on the D/H evolution in 
the local disk, based on the theoretical scheme that con- 
siders the wind enrichment and the decreasing SFR with a 
rate of ~ -0.8 Mq Gyr -1 , two phenomena that are implied 
to have occurred in the last several Gyr from the observed 
feature of local disk stars. The essence of its theory is that 
all heavy elements are not locally produced but some frac- 
tion comes form the outside, thus the models whose mission 
is the construction of the metal-rich system accept the ob- 
served low activity of star formation in the recent Gyrs, and 
end up with the ISM rich in D. 

We insist that the local star formation does not suffice 
to produce all elements in-situ, and thus D does not destroy 
so much as expected by the chemical features of the local 
disk. Si milar claim may be se e n in the mechanism of radia l 
mixing (|Roskar et al.l l2008bl ; ISchoenrich fe Binnevl [20091) , 
which brings metal-rich stars from th e outside (i.e., the in- 
ner d isk) instead of heavy elements. ISchoenrich fe Binnevl 
(2009) show the star formation history for the local disk 
which exponentially decreases 1 Gyr after the beginning of 
star form ation, as implied by the new reduction of Hippar- 
cos data (|Aumer fc Binnevl [2009). Their model succeeds in 
a high-Z stellar abundance system as an outcome of radial 
migration of stars while it will surely predict a low D today. 

We strongly predict an ongoing infall onto the Galaxy 



disk. Our claim is that it operates to reduce the gas abun- 
dances in the disk. Indeed, today's observed gas abundances 
are not metal-rich at any places of the disk. The metallici- 
ties of Cepheids at the i nner disk (R ^ 4kpc) are [Fe/H]~ 
+0.2 - 0.3 at the most (jAndrievskv et al.|[20o3 ). Moreover, 
a surprisingly rsult is brought by supergiants in two young 
clusters in the innermost disk at J?=3-4 kpc, whose metallic - 
ities are found to be [Fe/H]~ -0.1- -0.2 (|Davies et alj|2009h . 
An infall must rain down on today's disk. 
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